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Local unitary invariance and the notion of negativity fonts are used as the principle tools to 
construct four qubit invariants of degree 8, 12, and 24. A degree 8 polynomial invariant that is 
non-zero on pure four qubit states with four-body quantum correlations and zero on all other states, 
is identified. Classification of four qubit states into seven major classes, using criterion based on the 
nature of correlations, is discussed. 



To detect and quantify entanglement of composite quantum systems is a challenge taken up with great zeal by 
£\j ■ theorists and experimentalists alike. On the way, from the elegant bipartite separability criterion of Peres [l| up 
to classification schemes for four qubit states [2h11| . several useful ent ang lement measures and invariants have been 
found [12h31|. Two qubit entanglement is quantified by concurrence [32j, which for a pure state is equal to global 
negativity [12j . Entanglement of a three qubit state due to three-body quantum correlations is quantified by three 
tangle [lj] . For the most general three qubit state, the difference of squared global negativity and three tangle is a 
measure of two qubit correlations and satisfies CKW inequality [bij . A natural question is, which polynomial function 
of the coefficients quantifies entanglement due to four-body correlations? Can we write an invariant analogous to 
global negativity for two qubits and three tangle for three qubits to quantify four-body correlations? 



**i* > ', Invariant theory describes invariant properties of homogenous polynomials under general linear transformations. If 



we write a qubit state in multilinear form, we can find the set of invariants of the form in terms of state coefficients 
0'i 1 i 2 ...i N by using standard methods, as has been done in 0, [H, [l9|]. One may then investigate the properties of all 
invariants in the set. Our general aim, however, is to construct those polynomial invariants that quantify entanglement 
O" 1 due to if— body correlations in an N— qubit (N > K) pure state. This is done by constructing N— qubit invariants 
from multivariate forms with {K — 1) —qubit invariants as coefficients instead of cti 1 i 2 ...ij V - In particular, the invariant 
' that quantifies entanglement due to TV— body correlations is obtained from a biform having as coefficients the N — 1 
, qubit invariants. The term N— body correlations refers, strictly, to correlations of the type present in an A^— qubit 
GHZ state. The advantage of our approach [13, [H| is twofold. Firstly, we can choose to construct invariants that 
contain information about entanglement of a part of the system. Secondly, since the form of A^— qubit invariants is 
directly linked to the underlying structure of the composite system state, it can throw light on the suitability of a given 
state for a specific information processing task. Local unitary invariance and the notion of negativity fonts are used 
as the principle tools to identify AT— qubit invariants in an A^— qubit state. Negativity fonts are the elementary units 
of entanglement in a quantum superposition state. Determinants of negativity fonts are linked to matrices obtained 
from state operator through selective partial transposition (33l. l3~il| . In this article, we obtain analytical expressions for 
. . 1 polynomial invariants of degree 8, 12, and 24 for N = 4 states. One of the four qubit invariants is found to be non zero 
on states with four-body quantum correlations and zero on separable states as well as on states with entanglement 
due to two and three body correlations. It is analogous to three qubit invariant used to define three tangle [13| , and 
can likewise be used to construct an entanglement monotone to quantify four-body correlations. 

To obtain four qubit invariants that quantify four qubit quantum correlations, we follow a sequence of steps as 
given below: 

1. Identify two qubit invariants for a given pair in a three qubit state. 

2. Obtain a quadratic equation with two qubit invariants for a given pair of qubits as coefficients. Discriminant of 
the form is the three qubit invariant written in terms of two qubit invariants. 

3. Identify two qubit invariants in a four qubit state. Select three qubits and write three qubit invariants for these 
in a four qubit state. We identify five invariants, including two invariants analogous to ones known for a three qubit 
state. 

4. A local unitary on fourth qubit yields transformation equations for three qubit invariants. Proper unitaries 
can reduce the number of three qubit invariants in the set to four. The process of finding such local unitaries yields 
a quartic equation from which four qubit invariants are obtained. Since the invariants in a larger Hilbert space 
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are written in terms of relevant invariants in subspaces, it is possible to differentiate the invariants that quantify 
three-body quantum correlations from those that quantify four-body quantum correlations. 

In principle the process can be carried on to higher number of qubits. Polynomial invariants introduced by Luque 
and Thibon got geometrical meaning in the work of Levay [28j]. We point out the relation of our four qubit 
invariants with invariants in [l9| and (28[. 

Polynomial invariants that identify the nature of correlations in a state are useful to apply classification criteria 
proposed in [TTI j to four qubit states. Two multi qubit pure states are equivalent under stochastic local operations and 
classical communication (SLOCC) [HI if one can be obtained from the other with some probability using only local 
operations and classical communication amongst different parties. SLOCC equivalence is the central point in four 
qubit state classification into nine families in 2] . Borsten et al. [H have invoked the black- hole-qubit correspondence 
to derive the classification of four-qubit entanglement. However, it has been found that the number of four qubit 
SLOCC entanglement classes is much larger Q. The main result of Lamata et al. Q is that each of the eight genuine 
in-equivalent entanglement classes contains a continuous range of strictly non equivalent states, although with similar 
structure. O. Viehman et al. [13] select a set of generators for the SL(2, C)<g>4 -invariant polynomials or tangles 
and classify the eight families of ref. [3] using tangle patterns. In our classification scheme using correlation based 
criterion , multipartite states within the same class have same type of correlations but may have different number 
and type of negativity fonts in canonical state (all the states may not be SLOCC equivalent). In section IV, we 
calculate the relevant invariants for SLOCC families 0] and re-classify the states on the basis of number and nature 
of negativity fonts with non-zero determinants. The polynomial invariants used to classify the states in our scheme 
quantify correlations generated by distinct interaction types. For example, in a state with four-body correlations, 
each qubit must have interacted with all the others during state generation. Intutively, this information should 
be extremely useful to quantum state engineering. Negativity font analysis can be a helpful tool to optimize the 
subsystem interactions to tailor the invariant dynamics for a specific quantum information processing task. A minor 
point that will be discussed relates to the controversy regarding the family L a b 4 which is pointed out in ref. [29| to be 
a subclass of L a b C3 with (a = c), while in Q it has been shown that L Q {, 4 and \j a bc 3 belong to distinct SLOCC classes. 



I. NEGATIVITY FONTS AND TWO QUBIT INVARIANTS 

In this section, we briefly review the concepts of glob al partial pranspose [l|, global negativity [HI, [l8j], if— way 
partial transpose [35| and if— way negativity fonts [221 l23j |. We also identify those two qubit invariants which determine 
the entanglement of a pair of qubits in a three qubit state. 

A general N-qubit pure state may be written as 

|¥ Al ' A *- A ") = a ni2-m \hi2...i N ) , (1) 

where |ii«2---*jv) are the basis vectors spanning 2 N dimensional Hilbert space, and A p is the location of qubit p. The 
coefficients cii 1 i 2 ...i N are complex numbers. The local basis states of a single qubit are labelled by i m = and I, where 
rn = 1, N. The global partial transpose of an N qubit state p = \\fr A i, A 2,---A N ^ ^a 1 ,a 2 ,...a n | w ^ reS p ec t to qubit 
at location p is constructed from the matrix elements of p through 

(iii 2 — iwl P% Ap \h32—3N) = {iii2---ip-ijpip+i---iN\p\jii2---i P -iipj P +i--iN) ■ (2) 
To construct a if— way partial transpose [35| . every matrix element (iiia.-.ijvl p |jij2---iAr) is labelled by a number 

N 

if = (1 — ^im,im)> wnere &i m ,j m — 1 f° r **m — jm, arid <5j mJm — for i m ^ j m . Matrix elements of state operator 

m— 1 

with a given K represent if —way coherences present in the state. Local and non-local operations on a quantum 
superposition transform if —way coherences to if ± 1 way coherences. The if —way partial transpose of p with respect 
to subsystem p for if > 2 is obtained by selective transposition such that 

(hi2--iN\ p T K F \hh—3N) = (iii2—ip-ijpip+i—iN\p\ji32—j P -iipjp+i—jN), 

N 

if E( 1 -^.^) = X ' and S *pj, =0 ( 3 ) 

rn— 1 
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and 



{iii 2 ...i N \ p K p \j1j2-jN) = (iii2—iN\p\hh—3N) 



N 



(4) 



m— 1 



while 



(«i«2-iiv| P% P |jij2— jjv) = {hi2-ip-ijpip+i-iN\p\jij2-jp-iipj P +i--jN) 



N 



if E( 1 -^.^) = lor2 ' and *,j,=° 



(5) 



m — 1 



and 



(iii 2 ...i7v|p 2 p \hh-3N 



{lH2-lN\P\n32~3N) , 

if ^(l-<S Wm )^lor2. 



m— 1 

One can verify that global partial transpose may be expanded as 

AT 



Pg Ap = J2Pk p ~(N-2)p. 

K=2 



(6) 



(7) 



Negativity of p Tp , defined as N J 



P 



L — l) ,where \\p\\ 1 is the trace norm of p, arises due to all possible 



negativity fonts present in p Tp . Since p is a positive operators, global negativity depends on the negativity of if— way 
partially transposed operators with if > 2. 

To understand the concept of a negativity font in the context of an N— qubit system, consider the state 



A X A 2 ...A^ 
K 



\ = a ili2 ,,,i N \hi2--lN) + CLii+l,i 2 ...i N \h + hh—lN) 
+ a hh-3N \hh-3N) + a ji + l,j 2 ..-3N \h + 1,J2--JJV) 



N 



with if = (1 — &i m jm) an£ f — 0- The state 

m—l 



^AiA 2 ...A t 



N — if qubit product state. Let o^ 1 be the if— way partial transpose of ok 

sta 

and f^ 1 with negative eigenvalue given by 



is the product of a if —qubit GHZ-like state with 

with 



^MA 2 ...A t 



^MA 2 ...A t 



respect to qubit A±. If p is a pure state given by p = j^^-i^— Aw^ (^,AiA 2 ...a n | ^ ^gjj is a 4 x 4 sub-matrix of 



dot 



a i\i 2 ...i N a ji+l,j 2 ---ji' 
a il + lt 2 .. in a 3l32 --3N 



The matrix 



is referred to as a if —way negativity font [22|, [23|. A symbol used to represent 



a i\i 2 ...i N a ji+lj2---jf 
a ii+l,i 2 ...i N a 3\h---3N 

a negativity font, must identify the qubits that appear in if qubit GHZ-like state. Therefore we split the set of TV 
qubits with their locations and local basis indices given by, T — {(AiJ^ (A 2 ) i2 ... (Ajv) ijv } , into two subsets, with 
Si t containing qubits with local basis indices satisfying Si j — (i m ^ jm), and S2T having qubits for which 

Si mjm = 1 ( 

i m — jm)- To simplify the notation, we represent by s\ t T, the sequence of local basis indices for qubits in 
Si t- A specific negativity font is therefore represented by 



.»1»2 — «JV 



a i\i 2 ...i N a ji + l,j 2 ---jp 
a ii+l,i 2 ...i N a 3lh---3N 



A nonzero determinant Dg^ T T = det 

("sTt "' tN ) ensures that a K 1 is negative. A measurement on the state of a qubit 
with index in S±^t reduces ax to a separable state, whereas, measuring the state of a qubit in S^t does not change 
the negativity of <7^- Al • Elementary negativity fonts that quantify the negativity of p TAp for p ^ 1 are defined in 
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an analogous fashion. The determinant of a K— way negativity font detects K— body quantum correlations in an N 
qubit state. For even K, proper combinations of determinants of K— way negativity fonts are found to be invariant 
under the action of local operations on K qubits [23[ ■ 

For a two qubit state negative eigenvalue of partial transpose is the invariant that distinguishes between the 
separable and entangled states. Global negativity of | \J/j4i^2 ^ _ ^a ili2 is given by I AlA2 — |<Zooflii ~ aoi a io|> 

which is invariant under U Al ® U A2 . Here U Ai is a local unitary operator that acts on qubit Aj. The subscript on 



tAxAi 



refers to two-body correlations. A two qubit state therefore has a single negativity font v m = 



a general three qubit state, 



a oo a oi 
<2io an 



In 



A 1 A 2 A 3 \ _ 



> = E 



a ili 2 i 3 1*1*2*3/ 



the number of two-qubit invariants, for a selected pair of qubits, is three. For the pair A]A 2 j f° r example, these are 
determinants of 2— way negativity fonts defined as 



D 



and the difference (D mo - D° 



oo 

'(As 
( D 000 _| 



det 



a 00i 3 O01i 3 
Ol0i 3 0.1113 



, h = 0, 1, 



(8) 



D {> 



where 



D 0l2 ° = det 



a oi 2 o aoi 2 +i,i 

Oli 2 a l,i 2 +l,l 



12 = 0,1, 



(9) 



is determinant of a three-way negativity font. 



II. THREE-BODY CORRELATIONS AND THREE QUBIT INVARIANTS 



Our method was applied in ref. [22j , to construct three-tangle [1J| and a degree two four qubit invariant which is a 
function of determinants of 4-way negativity fonts. To clarify the process, we review the three qubit case and show 
that by using three qubit invariants one may classify three qubit entangled states into states with i) three and two 
body correlations, ii) states with only three body correlations and iii) a set of states with only two body correlations. 
Class (i) states are the most general states. Class (ii) states with GHZ type entanglement have the form 



AxA^Az 



) = ahi 2 i 3 1*1*2*3) + aj 1+M2+M3+ i \ii + + 1,«3 + 1) 



and Class (iii) contains W-like entangled states and biseparable states of three qubits. First of all, we write down the 
transformation equation for two qubit invariant D 9^ n to obtain the invariant which quantifies three-body correlations. 
The form of this invariant is later used to identify three qubit invariants in four qubit states. In the absence of three- 
body correlations, modified transformation equations yield three qubit invariants that quantify two body correlations 
in a three qubit state. 

"l -x* " 
x 1 



Under a local unitary \J J 



noo 

U (A 3 



transforms as 



(As), 



1 



(^ 3)i + (x) 2 ^ 3)o+ x(^ + ^ 01 )), 



(10) 



such that 



n oo 



U (A 3 



D 



000 



D 



001 



(11) 



is a three qubit invariant. If the pair of qubits A\ A2 is entangled then N^ A2 7^ 0. We can verify that global negativity 
of i^Q 1 is given by 



(^) 2 =4(<^) 2 +4(<^) 



where 



(12) 



D 



(A 2 ) ± 



D 



(A2) c 



/ J)000 _ J)001 \ 



(13) 
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The discriminant of ( ^(^3) ) = ^ ! y^l^s three qubit invariant 

T A 1 A 2 A 3 _ /£)' 



'(3,4) 



000 , r>ooi\ 2 a n°o r>oo 



(14) 



(15) 



which is a polynomial invariant of degree four in coefficients a ili2i3 . The subscript in J AlA2A3 refers to three-body 
correlations of the type present in a three qubit GHZ state. The terms D 00Q — D 010 , -D?^ 3 ) : an d -D?a 3 ) van ish on a 
product state of qubits A% and A 2 . On the state 

\^ AlM ) = Z ili2i3 *iiii l*i* 2 ) (bo |0) + 6x |1» ; (i™ = 0, 1) 

T AiA 2 A 

is a well known entanglement monotone. 

H 

1 



with D QQ £ 0, we have = (b fD 00 , = (h) 2 L> 00 , and L> 000 = D 0Q1 = b b x D m as such I A ^ A3 = 0. 

Moduh 
r 3 = 4 



Modulus of I 3 1 2 3 , quantifies the entanglement of qubits A1A2A3 due to three body correlations. Three tangle [14j . 



t a 1 a 2 a 3 
3 

For a general three qubit state with j AlA2A3 = 0, determinants of two-way fonts transform as 



( 



n oo 



1 + 



2 { x *\j~^X) 1 \I d Tm) ( ) ' 



noo 



therefore 



N 



A X A 2 



n oo 



n oo 



(16) 



is a three qubit invariant. In other words if l AlA2A ' 3 — Q then N AlA2 quantifies two body correlations of the pair 



A\A2. One can verify that 
product invariants 



D oo 



\A T , 



3 / 



(m = 1,2,3), are three qubit invariants in this case. The sum of 



D 



(Ai) 



+ 



D' 



(^)i 



)( 



D 



(Aj) 



D 



A\A 3 »,AiA : 



3 f N A,A 3N MA 3 + N MA 3N MA 2 + N MA 3N Ai 



(17) 



detects W-like tripartite entanglement. It is zero on bi-separable states for which only one of the three N 



AiAi 



noo 

n (A m ) 



+ 



noo 

U (A m ) 1 



(j ^ j ^ to) is non zero and one on a three qubit W-state. Major classes of three qubits states 

2 



are uniquely defined by values of polynomial invariants 4 



a ndI AlA * A \ 



III. FOUR-BODY CORRELATIONS AND FOUR-QUBIT INVARIANTS 

Four qubit states live in the Hilbert space C 2 (g) C 2 <E) C 2 ® C 2 with a distinct subspace for each set of three 
qubits. If there were no four body correlations, three qubit invariants ^j AzA ^ Ak ^j (j t = 0,1), may determine 

1 

the entanglement of a four qubit state. In general, additional three qubit invariants that depend also on four-way 
negativity fonts exist. For a selected set of three qubits, three qubit invariants constitute a five dimensional space 
and are easily found by the action of a local unitary on the fourth qubit. To write down transformation equations for 
three qubit invariants, first of all, we identify two qubit invariants. 
In the most general four qubit state 

\* AMsA *) = £ ili2i3i4 <W3H \ki2ki4) ; (»m = 0,l), (18) 



6 



when state of qubit Ai is transposed, we are looking at entanglement of qubit A\ with rest of the system. Qubit 
Ai may have pairwise entanglement with qubits A 2l A 3 , or A4. For a given pair, there are four two-way two qubit 
invariants (the remaining pair of qubits being in state |00), |10), |01) or |11)). For example, the determinants of 
two-way negativity fonts for the pair A\A 2 , written as 



a 00i 3 i 4 a 0l2 3 2 4 
&IO2324 ^112324 



, (i 3 ,i 4 = 0,1). 



(19) 



are invariant with respect to unitaries on qubits A\ and A 2 . Three-way coherences generate two qubit invariants 
D °(Z)- - D °(A 4 ). » (*4 = 0,1), and D°$ 3) - D° ( f 3) (i 3 = 0,1), for the pair A X A 2 . Here determinants of three-way 
fonts for {A!^^} and {Ai^^}, respectively, are defined as 



det 



«li 2 0i4 a \i 2 + \.\i& 



and 



(i 2 ,i 4 = 0, 1) , 



(20) 



det 



a 0i 2 i30 a 0i 2 + li3l 

Q*li 2 i 3 Q ^li 2 + l23l 



If four-way negativity fonts are present, then additional A\A 2 invariants, D 0000 
considered. Determinants of four-way negativity fonts are given by 



(i 2 ,i 3 = 0,1). 

D Q1QQ and D 0QQ1 



(21) 

D 0101 , are to be 



D 



Degree two four qubit invariant 



0i 2 0*4 



det 



O0i 2 0u a 0,i 2 + l,M 4 + l 
Oli 2 0i4 a l,i 2 + l,l : i4 + l 



(*2, *4 = 0, 1) . 



0000 



D' 



0011 



D 



0010 



D 



0001\ 



(22) 



(23) 



obtained in [22J is the same as invariant H of degree two in ref. [19J. The entanglement monotone, T4 = 4 1 J4 1 , was 
called four-tangle in analogy with three tangle [HI]. In [23| our method was successfully applied to derive degree two 
N— qubit invariants for even N and degree four invariants for odd N in terms of determinants of negativity fonts. 
It was also shown that one may use the method to construct N— qubit invariants to detect M— qubit correlations 



(M < iV) in an N— qubit state. As an example, we reported degree four invariants J, 



ref [23J and found that 



(/ 4 ) 2 = i(jf^)+ji^ + jf^). 



(A 1 A 2 ) r(AiA 3 ) 
) J 4 



and J, 



(AiA 4 ) 



Presently, we focus on the set AiA 2 A% of three qubits in state |vj/- 4 i j4 2- 4 3-4 4 ^ ([18]) ) viewed as 



A1A2A3A 



where 



tf,A 1 A 2 A 3 



4 ) 



**£*)| ) 



, T/ AiA 2 A 3 



(24) 



i 2 i 3 a i\i 2 i 3 ii 



\ki2i3i4)} («4 = 0, 1). 



Three qubit invariants 



A 2 A; 



/(A 4 ), V 



£)000 

(^4)i 



D 



001 



4n oo n oo 

iU (A 3 ) (Ai) z LJ {A 3 ) l {A i ) z 1 



i 4 = 0,1, 



(25) 



quantify GHZ state like three-way correlations in three qubit subspace C 2 <E) C 2 <E) C 2 . Continuing the search for 
a four qubit invariant that detects four qubit correlations, we examine the transformation of three qubit invariant 



fjAiAzAs^ 



(Ai), 



under 



\A+M 2 



1 -y* 
y 1 



A 2 A, 



The resulting transformation equation is 
1 



(A 4 



-6y 2 T^ A ^ 



1 (1 + M 2 )' 

>AiA 2 . 



A y P MA 2 A 3 



ifjA lA2 A 3 \ -4y 3 P^ 

{it- 



iA 2 A : 



(26) 



7 



where 



rpA 1 A 2 A 3 _ }_ /^OOOO _|_ £)0011 _|_ £,0010 , jjOOOIN 2 



6 
2 



and 



-(D ma +D oal ) (D mo +D 001 ) 

2 f n 00 n 00 -t- n 00 r> 00 ^ w\ 

3 \ U (A 3 ) (A i ) 1J (A 3 ) l (A i ) 1 + L '{A 3 ) a (A i ) 1 L, {A 3 ) 1 (A i ) ) > (/'J 



0A1A2A3 _ 1 / ,-,000 1 r)001 \ / 7-)0000 , rjOOll , rjOOlO , tjOOOI^i 



- 0(1)^), 0(A 3)o + ^ 3)o ) + ^l) (A 4)i4 + ^ 3)l )) ■ (28) 

Discriminant of a quartic equation, y 4 a — 4&y 3 + 6y 2 c — 4dy + / = 0, in variable y is A = S* 3 — 27T 2 where 
5 = 3c 2 — 4bd + af, and T = acf — ad 2 — b 2 f + 2bcd— c 3 (cubic invariant ), are polynomial invariants. When a selected 

\] Ai results in \\[l AlA2Aa ^ ^ ^ = (Eq. ([26]) )) the associated polynomial invariant is 



jA 1 A 2 A 3 A 4 _ 3 /2iAiA 2 A 3 \ 2 _ ^pAiA 2 A 3 pAiA 2 A 3 
(4,8) \ A 4 J (A 4 ) (Ai)i 



V 3 V 3 /(M 



(A 4 ) V / (Ai)i 

which is a four qubit invariant of degree eight expressed in terms of three qubit invariants for A1A2A3. In order 
to distinguish between degree 2 invariant I 4 and the new invariant, degree of the invariant has been added to the 
subscript. By construction, the four qubit invariant I^ 1 ^ 2 ^ 3 ^ is a combination of three qubit (A1A2A3) invariants. 
It is easily verified that on a state which is a product of ^-^i^a-As^ = £\ iai a^i.,^ 1*1*2*3) with lAiA 2 A 3 _^ &nd 
H Ai = do |0) + di |1) , we obtain 



{ jA 1 A 2 A 3 \ _ ( jA!A 2 A 3 \ 



(A 4 )! 



_ ^4.1^2^3 

l Ai 



p, 



p, 



(30) 



leading to i^ 1 ^ 3 ^ = 0. Likewise, I^ 2 ' 4 ' 3 ' 44 vanishes on product state \$ AlA2 ) ^A 3 A A ^ w here |^ Al/l2 ) = 



E ili2 aiii2 |»i*2> and | = £^6^ 1*3^4 

entangled states with only three and two-body correlations, as seen in section IV. 
The cubic invariant associated with Eq. (|2l)l) is 



Besides that I^g) = on a four qubit W-like state, and all 



J 



AiA 2 A 3 Ai = det 



A]_A 2 A 3 



P, 



A X A 2 A 3 



r>A\A 2 A 3 r r A 1 A 2 A 3 
r (A i ) 1 1 A l 



pA±A 2 A 3 
-A 4 



P, 



Ai_A 2 A 3 



rpA 1 A 2 A 3 

1 A i 

pA 1 A 2 A 3 
^4)o, 
A1A2A3 



(^4) 



(31) 



while the discriminant reads as 



A = 



(tAi 



A 2 A 3 Al 



27 (J 



A\A 2 A 3 An\ 



(32) 



Since there are four ways in which a given set of three qubits may be selected, A can be expressed in terms of different 
sets of three qubit invariants. In addition (Eq. ([2"B)) ) also leads to 



2 



+6 



Ti 



+ 4 



r>A 1 A 2 A 3 



{A i ) 1 



P, 



(33) 



8 



which is a four qubit invariant analogous to 

( N A 3 A2 ) (Eq- CU!) for three qubit states. In general, one can construct 

an invariant jSf AzA]Ak (j ^ j ^ f~ ^ /) for a selected three qubit subsystem AiAjAk of four qubit state. In analogy 
with global negativity, one may define a four qubit invariant of degree four, 

(iV ( t 4) ) 2 = 16 (Nil A ^) 2 + 16 (N A l A ^) 2 + 16 (N A l A ^)\ (34) 

which detects bipartite entanglement of qubit A\ with subsystem A2A3A4 due to three and four body quantum 
correlations. If I^g^ 3 ^ = 0, but at least two of the N A * A * Ak are finite, then 4-partite entanglement can be due to 
three and two body correlations. In this case the invariant that detects entanglement may be defined as 

N A ) A2MM = \QN A l MM N A l MM + 16 (N A f 3M + N Al 2 MA ^ N A l MM 

+ 16 (n a * mm + N^ AsAi + N A l MM ) N£ AaA * . (35) 

nOO 



On the other hand, if we have a state on which all N A ' 3 k are zero, then the quantities I J s L p J s L q = X) 
(p ^ q ^ r 7^ s = 1 to 4), turn out to be four qubit invariants. A different class of entangled states is obtained if only 

A - A ■ Ai A A 

I A % 3 ' is non zero along with a finite I £4 

'(4,8) - Jv (4,8) 



one of the AT? 3 is non zero along with a finite ■ In section II we noted that I A ^ As (Eq. (JYTJ)) detects W 

like entanglement of three qubits A 1 A 2 A 3 . Likewise, when I A ^ AsA * = N AlA2A ' iAi = 0, the invariant 



TA1A2A3A4 _ 3 / T AiA a A 3 \ I jAiAi , J-A2A4 , T A 3 A 
J (2,6) - 2 V 2 / V AlM + 1 A 1 A 

. 3 / T 4 , 4 „ 4 , \ ( T 4 „ A „ T 4,4,\ . 3 



(t a ^ a A (i 



2 

detects W-like four qubit entanglement. Here 

^ApAtAr^j _ 3/^y^;, (p ^ g ^ r ^ ,s = 1 to 4), is the invariant 

that detects W-like entanglement of qubits A p A q A r in a four qubit state. 

In ref. [28| four qubit invariants have been obtained in terms of coefficients having geometrical significance. A 
comparison of Eq. (56) of ref. [28] with our Eq. (|26l) . indicates that their set of invariants (ii, I2, I3, I 4) may be 
expressed in terms of our three qubit invariants, though they are not exactly the same. A method equivalent to 
method of Schlafli (36] has been used to arrive at Eq. (22) 

R(t) = cc4 + 4cit^i + 6c2tl4 + 4c 3 t *? + C4.4 

by Luque and Thibon (l9| . Then higher degree invariants are expressed in terms of Cj coefficients and computer 
algebra relates these to basic four qubit invariants. Since for to = 1, expression for R(t) has the same form as 
Eq. (|26[) . a direct correspondence can be established between Ci coefficients and our three qubit invariants. Such a 
comparison establishes a neat connection of our invariants with projective geometry approach and classical invariant 
theory concepts. 



IV. INVARIANTS AND CLASSIFICATION OF FOUR-QUBIT STATES 

Decomposition of global partial transpose p^ Ap of four qubit state ^^^'^Ai^ w ith respect to qubit A p in terms 
of K— way partially transposed operators (Eq. (J7J) reads as 

f G Ap = E P T k P ~ 2£ (37) 

K=2 

When a state has only K— way coherences, we have p^ v = 'p^ v , for a selected set of K qubits. For a given qubit, 
the number of if —way negativity fonts in a if —way partially transposed matrix varies from to 4. Local unitary 
operations can be used to annililate the negativity fonts that is obtain a state for which determinants of selected 
negativity fonts are zero. The process leads to canonical state which is a state written in terms of minimum number 



of local basis product states [37|. In ref. we proposed a classification scheme in which an entanglement class 

is characterized by the minimal set of if —way (2 < K < 4) partially transposed matrices present in the expansion 
of global partial transpose of the canonical state. Seven possible ways in which the global partial transpose (GPT) 



9 



TABLE I: Decomposition of (pc) G Ap , invariants I^^ 2A3Ai , D A * A2A3 , A, and iVx-way (K =2,3,4) in canonical state, for seven 
classes of four qubit entangled states 



Class 


Decomposition of (p c ) G Ap 


rA 1 A 2 A 3 A i 
J (4,8) 


D MA 2 A 3 


A 


N2 — may 


■^V 3 — way 


— way 


I 


J2 {Pc) T k V - 2 Pc 

if =2 









> 1 


> l 


> l 


TT 
11 


lPcJ 4 + (Pc) 3 - pc 


-/- n 
T u 




u 


A 
U 


^ 1 


1 


in 


(pc) 4 A " + (pc) 2 Ap ~Pc 









> 1 





> l 


IV 




¥=0 














l 


V 


(Pc)s AP + (Pc)^ - Pc 










> 1 


> 1 





VI 


(Pc)a Ap 













1 





VII 


(Pc)^ 











> 1 









of a four qubit canonical state may be decomposed correspond to seven major entanglement classes that is class 

I- (Pc)l Ap = £ (Pc) T K P ~ 2?c, II. (Pcf G Ap = {Pcf^ + (Pc)l Ap ~ Pc, HI- {Pcf G Ap = {Pcf^ + {Pcf 2 Ap ~ Pc , IV. 
K=2 

{Pcf G Ap = {Pc)l Av , V. (p c f G Ap = {p c )l Ap + {p c f 2 Ap -p c , VI. (p c f G Ap = {p c )l Ap , and VII. {pf G Ap = (p c f 2 Ap . Of 
these, six classes contain states with four-partite entanglement, while class VI with p^ Q Ap — (j^ Ap ^J has only three 
qubit entanglement. Each major class contains sub-classes depending on the number and type of negativity fonts 
in global partial transpose of the canonical state. Table Q] lists the decomposition of (/O c )g Ap , invariants I A ^ A } 2A3Ai , 
D AlA2A3 , A, and Nk -way {K =2,3,4) in canonical state, for different classes of four qubit entangled states. Here 

2 .... . 

is a measure of residual three-way correlations between qubits A1A2A3 and 



£ ) A 1 A 2 A 3 _ I AT A 1 A 2 A 3 



1 (4,8) 



^V/f-way (K =2,3,4) is the number of if —way negativity fonts in a state. 
A four qubit state with a single four way negativity font 

|* a6 ) = a (10000) + |1111)) + b (|1101) + |1110) + |0011» , 

is an example of class I states . Three qubit invariants for the state are lf} A2Aa = a 2 b 2 , Pf\ lA2Aa = ±a 3 6, I A , lA2A ' 3 = 
b 4 , Pf2f)l As = -\a 2 b 2 , and (T AlA * Aa ) Ai = ±(a 4 -2a6 3 ). Four qubit invariants are found to be I^ 1 ^ 3 ^ = 

i (a 4 + 4ab 3 ) 2 , D AlA2M ^ 0, and A / 0. A representative of class II states with p^ Ap = f Ap + p% Ap - p is, 
|* a ) = a (|0000) + 1 1 1 1 1>) + |1110). The state is SLOCC equivalent to GHZ state, however it deserves a distinct 
status since on removal of qubit A4 it has residual three way coherences. 

Invariants for class III states G abcd , L abC2 , i a2&2 and L a2 o 3eI with (p c )c Ap = (Pc)!"* 1 + (Pc)^ 1 - Pc are listed in 
Table [TT] All three way coherences are convertible to two way coherences as such three-way negativity fonts have zero 
determinants. Four qubit entanglement occurs due to four-way and two-way coherences. For all these states, the 
invariants P^l^ 2 ^ and P^^ 2 ^ are identically zero. In Table HU for states in family G a bcd (ref. 0), three qubit 
invariants used for the set A1A2A3 are 



A 2 M 



(^4)0 V / (A i ) 1 



B, 



where 



A = (a 2 - b 2 ) (d 2 -c 2 ),B = ±(a 2 - d 2 ) (b 2 - c 2 ) . (38) 



(A0 Q 



For states G ab00 and G a0cd , A = 0. For states L abC2 , with (T A ^ A ' A ") M = \{a 2 - c 2 ) (b 2 - c 2 ), (l AM ^ 

c (a 2 — fe 2 ), the value ( j^ lA2Aa J =0 results in A = 0. A comparison of states L abC2 with a — c and L ab3 shows 
that the states are not SLOCC equivalent [8j because the number of negativity fonts is not equal. However, since 
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TABLE II: Invariants for class III states G a bcd, L a bc 2 , L a2 b 2 an d L a2 o j with (pc) G p = (pc) 4 1 + (pc) 2 1 — pc- A and B in 
column II are as defined in Eq. ((38}. 



Invariant\Class 


Gabcd 








(K A2A3 ) 2 


| \A-2B\ 2 + 2\B\ 2 


11/2 „2\ /t2 2\ |2 

+ \c(a 2 -b 2 )\ 2 


i (a 2 -b 2 ) 4 




tA 1 A 2 A 3 A 4 
J (4,8) 


(±(A-2Bf+B*) 


±(a 2 -c 2 y(b 2 ~c 2 ) 2 


jr 2 {a 2 -b 2 Y 




jjA 1 A 2 A 3 
Ad 




c (a 2 - b 2 ) 2 








A 


7^0 












four qubit correlations are null (^(t^g) 2 " 43 " 44 = ®j lor £ahe 2 w hh a = cas well as L a f, 3 , these are subclasses of the same 
major class in correlation type based classification, partially supporting the result of [29]. 

The families of states L a b 3 and L ai (ref. [2]) have a similar composition of global partial transpose. The value of 
degree two invariant is I4 — 3a + 6 for L a (, 3 and I4 = 2a 2 for L ai indicating that four-way coherences are present. 
However, for the set of qubits only non zero three qubit invariant is (l AlA2Aa j = a ~ b for L a {, 3 

V / {A 4 ) 1 

and (l AlA2A3 j = —4a 2 for L ai . A finite I4 but zero I^ 1 ^ 2 ^^ indicates that the superposition contains 
a product of two qubit entangled states. Four partite entanglement may, in this case, be detected by products 

rA 1 A 2 A 3 A i 



The states in families L no ^ and L nn a . have 



„,,,.. aim i.„ | nave (n%1 



A 2 A 3 \ 2 _ 



TJ\\J\ 2 J\ 3 I\i 

J (4,8) 



. The states in L a2 b 2 and £ a2 o _ differ 



from each other in the number of two way negativity fonts with non-zero determinants. Only non zero three tangle 
for the states G a bba is T^ A2As .The states G a ooa and Gobbo, with (p c )q Ap — (/Oc)J Ap belong to class IV in classification 
scheme based on correlation type. For these states only non-zero three tangle is 3 , therefore I A ^ A2A:>Ai 7^ 0, 

A = and D AlMM = 0. 

A4 

The global partial transpose has composition, (p c ) G Al = {pc)^ 1 + (pc) 2 Al — Pc, for class V states Lo T and £o 5(B3 - 
In both cases I^ MM = 0, while the product (a^ 1 ^ 3 ) (a^ 3 iA2j44 ) ^ 0. Two states differ in the the number 
of two-way negativity fonts with non-zero determinants. Only non-zero invariant for Class VI state Lq _o _,with 

3© 1 3Q-) 1 

(P^lc? 1 = {Pc)^ Al is [I3 2 3 4 ) ■ The state has only three qubit entanglement. Class VII with (p c ) G p — {Pc)^ 

contains four qubit states with W-type entanglement represented by L a= of, 3= o and separable states with entangled 
qubit pairs, for example G aaaa . 

The polynomial invariant I'u^) 2A ' iAi is non-zero on states \^ a b), G a bcd, L a bc 2 , L a2 b 2 , L a2 o ~, G a oo a and Gobbo and 



\ Ta v 
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vanishes on states L ob3 , L ai , io 7@T , ^o 5ffi3 , L o mT o 3S) j and G aaaa . We define an entanglement monotone to quantify 



'(4,8) uuu-aciu un aucxuco \*at>/, ^abca, -^abc 2 , -^a 2 b 2 , -^a 2 U 3e 

four qubit correlations as 

J A 1 A 2 A 3 A i 



T(4,8) = 4 



197" 

(4,8) 



(39) 



which is one on states with maximal entanglement due to four-body correlations, finite on all states with entanglement 

T A t A 2 ' ' 
'(4,8) 



due to four-body correlations and zero otherwise. The subscript (4,8) is carried on from -/J 4 ^ 2 " 43 " 4 ". One can verify 



that on four qubit GHZ state 



as well as cluster states [38!, |3£ 



\GHZ) = -j= (10000) + 11111)) 



\Ci) = i (10000) + |1100) + |0011> - |1111» 



|C 2 ) = - (10000) + 10110) + |1001) - |1111» , 
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\c 3 



- (|0000) + |1010) + |0101) - |1111)) . 



T(4,8) 



1 and 



(< 



A 2 A- 



TA1A2A3A4 
J (4,8) 



So what is different in cluster states? We recall the invariants J AiA i 



from [23j, the invariants that detect entanglement of a selected pair, AjAj, of qubits in a four qubit state. For a GHZ 
state J AiA i = ij for (i ^ j) = 1 to 4, while for a cluster state all J AiA i 23], do not have the same value. In canonical 
form, GHZ has a single four-way negativity font, while a cluster state has two four-way negativity fonts besides also 
having two-way negativity fonts (state reduction does not destroy all the coherences). 

Another state proposed through a numerical search in ref. [40] to be a maximally entangled state is 



|$) = - (10000) + |1101» + -= (|1011) + 10011) + |0110) - |1110» . 
2 v 8 



However, on this state 



iA 1 A 2 A 3 



A 2 A 3 



(^4)0 



= I 



r A 1 A 2 A i 



(A 4 



1 

32' 



f pA 1 A 2 A 3 \ _ (pA t A 2 A 3 \ 

\ 3 )(A i ) a \ 3 ){A i ) 1 



0. 



therefore lfc£» AaA 1 - 



— ^g, and 774,8) =yf- On two excitation four qubit Dicke state 



\9d) = -7= (10011) + |1100) + 10101) + |1010) + |1001) + |0110)) , 
V6 

we have, 774,8) — § j while it is zero on four qubit W-state 

\W) = i (|0000) + 1 1100) + |1010) + 1 1001)) . 

Four tangle 74 also vanishes on W— like state of four qubits, however, it fails to vanish on product of two qubit 
entangled states. Contrary to T( 4 8 ), a non zero 74 does not ensure four-partite entanglement. On four qubit state 

\HS) = i (J0011) + 11100) + exp (jp) (11010) + |0101» 

+ ± eX p(^\ (11001) + 0110), (40) 

conjectured to have maximal entanglement in ref. [4l|, we have -D™ 3 ) <Ai) = ^(A 3 ) = I' an< ^ ^ or 4-way 

negativity fonts L> 0011 = \, D mm = i (l _ iy^) , and D 0010 = i (l + ^)). Therefore 



T 



^1^2,43 _ / jA\A 2 A 



(^4) 



_ (jA x A 2 A 3 



(A i ) 1 



= 0, 



V )( Ai ) n \ JiAi), 

leading to 774,8) = 0. However, the invariant I^ 1 ^ 2 ^^ = 1 (Eq. dM])) on \HS) and takes value || on four qubit |W) 
state. It reflects the fact that a measurement on the state of a qubit, in \HS) always leaves the three remaining qubits 
in a three qubit W-state, whereas a similar measurement on a | W) state yields a mixture of three qubit W-state with 
three qubits in a separable state. 

The choice I^ 1 ^ 2 ^^ to quantify four qubit correlations is also supported by the conclusions of 42], where for a 
selected set of four qubit states, generator S of ref. has been shown to have the same parameter dependence as 
optimized Bell type inequalities and a combination of global negativity and 2-qubit concurrences. 

To summarize, degree 8, 12 and 24 four qubit invariants, expressed in terms of three qubit invariants, have been 
obtained. One can continue the process to higher number of qubits. Commonly, multivariate forms in terms of state 
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coefficients ai t i 2 ,,,i N are used to obtain polynomial invariants for qubit systems. Our strategy is to write multivariate 
forms with relevant K— qubit invariants as coefficients. The advantage of our technique is that relevant invariants in a 
larger Hilbert space are easily related to invariants in sub spaces as such to the structure of the quantum state at hand. 
Construction of polynomial invariants for states other than the most general state is a great help in classification of 
states. Our method can be easily applied to determine the invariants for any given state. Entanglement monotone that 
quantifies four qubit correlations can be used to quantify correlations in pure and mixed (via convex roof extension) 
four qubit states. 

This work is financially supported by CNPq Brazil and FAEP UEL Brazil. 
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